Validation of daily 0.35 T diffusion-weighted MRI for MRI-guided glioblastoma radiotherapy.

BACKGROUND: MRI-Linac systems enable daily diffusion-weighed imaging (DWI) MRI scans for assessing glioblastoma tumor changes with radiotherapy treatment. PURPOSE: Our study assessed the image quality of echoplanar imaging (EPI)-DWI scans compared with turbo spin echo (TSE)-DWI scans at 0.35 Tesla (T) and compared the apparent diffusion coefficient (ADC) values and distortion of EPI-DWI on 0.35 T MRI-Linac compared to high-field diagnostic MRI scanners. METHODS: The calibrated National Institute of Standards and Technology (NIST)/Quantitative Imaging Biomarkers Alliance (QIBA) Diffusion Phantom was scanned on a 0.35 T MRI-Linac, and 1.5 T and 3 T MRI with EPI-DWI. Five patients were scanned on a 0.35 T MRI-Linac with a TSE-DWI sequence, and five other patients were scanned with EPI-DWI on a 0.35 T MRI-Linac and a 3 T MRI. The quality of images was compared between the TSE-DWI and EPI-DWI on the 0.35 T MRI-Linac assessing signal-to-noise ratios and presence of artifacts. EPI-DWI ADC values and distortion magnitude were measured and compared between 0.35 T MRI-Linac and high-field MRI for both phantom and patient studies. RESULTS: The average ADC differences between EPI-DWI acquired on the 0.35 T MRI-Linac, 1.5 T and 3 T MRI scanners and published references in the phantom study were 1.7%, 0.4% and 1.0%, respectively. Comparing the ADC values based on EPI-DWI in glioblastoma tumors, there was a 3.36% difference between 0.35 and 3 T measurements. Susceptibility-induced distortions in the EPI-DWI phantoms were 0.46 ± 1.51 mm for 0.35 MRI-Linac, 0.98 ± 0.51 mm for 1.5 T MRI and 1.14 ± 1.88 mm for 3 T MRI; for patients -0.47 ± 0.78 mm for 0.35 T and 1.73 ± 2.11 mm for 3 T MRIs. The mean deformable registration distortion for a phantom was 1.1 ± 0.22 mm, 3.5 ± 0.39 mm and 4.7 ± 0.37 mm for the 0.35 T MRI-Linac, 1.5 T MRI, and 3 T MRI scanners, respectively; for patients this distortion was -0.46 ± 0.57 mm for 0.35 T and 4.2 ± 0.41 mm for 3 T. EPI-DWI 0.35 T MRI-Linac images showed higher SNR and lack of artifacts compared with TSE-DWI, especially at higher b-values up to 1000 s/mm2. CONCLUSION: EPI-DWI on a 0.35 T MRI-Linac showed superior image quality compared with TSE-DWI, minor and less distortions than high-field diagnostic scanners, and comparable ADC values in phantoms and glioblastoma tumors. EPI-DWI should be investigated on the 0.35 T MRI-Linac for prediction of early response in patients with glioblastoma.

Performance of deep learning synthetic CTs for MR-only brain radiation therapy.

PURPOSE: To evaluate the dosimetric and image-guided radiation therapy (IGRT) performance of a novel generative adversarial network (GAN) generated synthetic CT (synCT) in the brain and compare its performance for clinical use including conventional brain radiotherapy, cranial stereotactic radiosurgery (SRS), planar, and volumetric IGRT. METHODS AND MATERIALS: SynCT images for 12 brain cancer patients (6 SRS, 6 conventional) were generated from T1-weighted postgadolinium magnetic resonance (MR) images by applying a GAN model with a residual network (ResNet) generator and a convolutional neural network (CNN) with 5 convolutional layers as the discriminator that classified input images as real or synthetic. Following rigid registration, clinical structures and treatment plans derived from simulation CT (simCT) images were transferred to synCTs. Dose was recalculated for 15 simCT/synCT plan pairs using fixed monitor units. Two-dimensional (2D) gamma analysis (2%/2 mm, 1%/1 mm) was performed to compare dose distributions at isocenter. Dose-volume histogram (DVH) metrics (D95% , D99% , D0.2cc, and D0.035cc ) were assessed for the targets and organ at risks (OARs). IGRT performance was evaluated via volumetric registration between cone beam CT (CBCT) to synCT/simCT and planar registration between KV images to synCT/simCT digital reconstructed radiographs (DRRs). RESULTS: Average gamma passing rates at 1%/1mm and 2%/2mm were 99.0 ± 1.5% and 99.9 ± 0.2%, respectively. Excellent agreement in DVH metrics was observed (mean difference ≤0.10 ± 0.04 Gy for targets, 0.13 ± 0.04 Gy for OARs). The population averaged mean difference in CBCT-synCT registrations were <0.2 mm and 0.1 degree different from simCT-based registrations. The mean difference between kV-synCT DRR and kV-simCT DRR registrations was <0.5 mm with no statistically significant differences observed (P > 0.05). An outlier with a large resection cavity exhibited the worst-case scenario. CONCLUSION: Brain GAN synCTs demonstrated excellent performance for dosimetric and IGRT endpoints, offering potential use in high precision brain cancer therapy.

Impairments of white matter tracts and connectivity alterations in five cognitive networks of patients with multiple sclerosis.

INTRODUCTION: MS is associated with structural and functional brain alterations leading to cognitive impairments across multiple domains including attention, memory, and speed of information processing. Here, we analyzed the white matter damage and topological organization of white matter tracts in specific brain regions responsible for cognition in MS. METHODS: Brain DTI, rs-fMRI, T1, T2, and T2-FLAIR were acquired for 22 MS subjects and 22 healthy controls. Automatic brain parcellation was performed on T1-weighted images. Skull-stripped T1-weighted intensity inverted images were co-registered to the b0 image. Diffusion-weighted images were processed to perform whole brain tractography. The rs-fMRI data were processed, and the connectivity matrixes were analyzed to identify significant differences in the network of nodes between the two groups using NBS analysis. In addition, diffusion entropy maps were produced from DTI data sets using in-house software. RESULTS: MS subjects exhibited significantly reduced mean FA and entropy in 38 and 34 regions, respectively, out of a total of 54 regions. The connectivity values in both structural and functional analyses were decreased in most regions of the default mode network and in four other cognitive networks in MS subjects compared to healthy controls. MS also induced significant reduction in the normalized hippocampus and corpus callosum volumes; the normalized hippocampus volume was significantly correlated with EDSS scores. CONCLUSION: MS subjects have significant white matter damage and reduction of FA and entropy in various brain regions involved in cognitive networks. Structural and functional connectivity within the default mode network and an additional four cognitive networks exhibited significant changes compared with healthy controls.

Rapid multicontrast brain imaging on a 0.35T MR-linac.

PURPOSE: Magnetic resonance-guided radiation therapy (MRgRT) has shown great promise for localization and real-time tumor monitoring. However, to date, quantitative imaging has been limited for low field MRgRT. This work benchmarks quantitative T1, R2*, and Proton Density (PD)mapping in a phantom on a 0.35 T MR-linac and implements a novel acquisition method, STrategically Acquired Gradient Echo (STAGE). To further validate STAGE in a clinical setting, a pilot study was undertaken in a cohort of brain tumor patients to elucidate opportunities for longitudinal functional imaging with an MR-linac in the brain. METHODS: STAGE (two triple-echo gradient echo (GRE) acquisitions) was optimized for a 0.35T low-field MR-linac. Simulations were performed to choose two flip angles to optimize signal-to-noise ratio (SNR) and T1-mapping precision. Tradeoffs between SNR, scan time, and spatial resolution for whole-brain coverage were evaluated in healthy volunteers. Data were inputted into a STAGE processing pipeline to yield four qualitative images (T1-weighted, enhanced T1-weighted, proton-density (PD) weighted, and simulated FLuid-Attenuated Inversion Recovery (sFLAIR)), and three quantitative datasets (T1, PD, and R2*). A benchmarking ISMRM/NIST phantom consisting of vials with variable NiCl2 and MnCl2 concentrations was scanned using variable flip angles (VFA) (2-60 degrees) and inversion recovery (IR) methods at 0.35 T. STAGE and VFA T1 values of vials were compared to IR T1 values. As measures of agreement with reference values and repeatability, relative error (RE) and coefficient of variability (CV) were calculated, respectively, for quantitative MR values within the phantom vials (spheres). To demonstrate feasibility, longitudinal STAGE data (pretreatment, weekly, and ~ 2 months post-treatment) were acquired in an IRB-approved pilot study of brain tumor cases via the generation of temporal and differential quantitative MRI maps. RESULTS: In the phantom, RE of measured VFA T1 and STAGE relative to IR reference values were 7.0 ± 2.5% and 9.5 ± 2.2% respectively. RE for the PD vials was 8.1 ± 6.8% and CV for phantom R2* measurements was 10.1 ± 9.9%. Simulations and volunteer experiments yielded final STAGE parameters of FA = 50°/10°, 1 × 1 × 3 mm3 resolution, TR = 40 ms, TE = 5/20/34 ms in 10 min (64 slices). In the pilot study of brain tumor patients, differential maps for R2* and T1 maps were sensitive to local tumor changes and appeared similar to 3 T follow-up MRI datasets. CONCLUSION: Quantitative T1, R2*, and PD mapping are promising at 0.35 T agreeing well with reference data. STAGE phantom data offer quantitative representations comparable to traditional methods in a fraction of the acquisition time. Initial feasibility of implementing STAGE at 0.35 T in a patient brain tumor cohort suggests that detectable changes can be observed over time. With confirmation in a larger cohort, results may be implemented to identify areas of recurrence and facilitate adaptive radiation therapy.

Large field of view distortion assessment in a low-field MR-linac.

PURPOSE: MR-guided radiation therapy (RT) offers unparalleled soft tissue contrast for localization and target tracking. However, MRI distortions may be detrimental to high precision RT. This work characterizes the gradient nonlinearity (GNL) and total distortions over the first year of clinical operation of a 0.35T MR-linac. METHODS: For GNL characterization, an in-house large field of view (FOV) phantom (60 × 42.5 × 55 cm3 , >6000 spherical landmarks) was configured and scanned at four timepoints with forward/reverse read polarities (Gradient Echo sequence, FA/TR/TE = 28°/30 ms/6 ms). GNL was measured in Anterior-Posterior (AP), Left-Right (LR), and Superior-Inferior (SI) frequency-encoding directions based on deviation of the auto-segmented landmark centroids between rigidly registered MR and CT images and assessed based on radial distance from magnet isocenter. Total distortion was assessed using a 30 × 30 cm2 grid phantom oriented along the cardinal axes over >1 year of operation. RESULTS: The scanner's spatial integrity within the first ~10 months was stable (maximum total distortion variation = 10/6/8%, maximum distortion = 1.41/0.99/1.56 mm in Axial/Coronal/Sagittal planes, respectively). GNL distortions measured during this time period <10 cm from isocenter were (-0.74, 0.45), (-0.67, 0.53), and (-0.86, 0.70) mm in AP/LR/SI directions. In the 10-20 cm range, <1.5% of the distortions exceeded 2 mm in the AP and LR axes while <4% of the distortions exceeded 2 mm for SI. After major repairs and magnet re-shim, detectable changes were observed in total and GNL distortions (20% reduction in AP and 36% increase in SI direction in the 20-25 cm range). Across all timepoints and axes, 38-53% of landmarks in the 20-25 cm range were displaced by >1 mm. CONCLUSIONS: GNL distortions were negligible within a 10 cm radius from isocenter. However, in the periphery, non-negligible distortions of up to ~7 mm were observed, which may necessitate GNL corrections for MR-IGRT for treatment sites distant from magnet isocenter.

Geometric and dosimetric impact of anatomical changes for MR-only radiation therapy for the prostate.

PURPOSE: With the move towards magnetic resonance imaging (MRI) as a primary treatment planning modality option for men with prostate cancer, it becomes critical to quantify the potential uncertainties introduced for MR-only planning. This work characterized geometric and dosimetric intra-fractional changes between the prostate, seminal vesicles (SVs), and organs at risk (OARs) in response to bladder filling conditions. MATERIALS AND METHODS: T2-weighted and mDixon sequences (3-4 time points/subject, at 1, 1.5 and 3.0 T with totally 34 evaluable time points) were acquired in nine subjects using a fixed bladder filling protocol (bladder void, 20 oz water consumed pre-imaging, 10 oz mid-session). Using mDixon images, Magnetic Resonance for Calculating Attenuation (MR-CAT) synthetic computed tomography (CT) images were generated by classifying voxels as muscle, adipose, spongy, and compact bone and by assignment of bulk Hounsfield Unit values. Organs including the prostate, SVs, bladder, and rectum were delineated on the T2 images at each time point by one physician. The displacement of the prostate and SVs was assessed based on the shift of the center of mass of the delineated organs from the reference state (fullest bladder). Changes in dose plans at different bladder states were assessed based on volumetric modulated arc radiotherapy (VMAT) plans generated for the reference state. RESULTS: Bladder volume reduction of 70 ± 14% from the final to initial time point (relative to the final volume) was observed in the subject population. In the empty bladder condition, the dose delivered to 95% of the planning target volume (PTV) (D95%) reduced significantly for all cases (11.53 ± 6.00%) likely due to anterior shifts of prostate/SVs relative to full bladder conditions. D15% to the bladder increased consistently in all subjects (42.27 ± 40.52%). Changes in D15% to the rectum were patient-specific, ranging from -23.93% to 22.28% (-0.76 ± 15.30%). CONCLUSIONS: Variations in the bladder and rectal volume can significantly dislocate the prostate and OARs, which can negatively impact the dose delivered to these organs. This warrants proper preparation of patients during treatment and imaging sessions, especially when imaging required longer scan times such as MR protocols.

Generating synthetic CTs from magnetic resonance images using generative adversarial networks.

PURPOSE: While MR-only treatment planning using synthetic CTs (synCTs) offers potential for streamlining clinical workflow, a need exists for an efficient and automated synCT generation in the brain to facilitate near real-time MR-only planning. This work describes a novel method for generating brain synCTs based on generative adversarial networks (GANs), a deep learning model that trains two competing networks simultaneously, and compares it to a deep convolutional neural network (CNN). METHODS: Post-Gadolinium T1-Weighted and CT-SIM images from fifteen brain cancer patients were retrospectively analyzed. The GAN model was developed to generate synCTs using T1-weighted MRI images as the input using a residual network (ResNet) as the generator. The discriminator is a CNN with five convolutional layers that classified the input image as real or synthetic. Fivefold cross-validation was performed to validate our model. GAN performance was compared to CNN based on mean absolute error (MAE), structural similarity index (SSIM), and peak signal-to-noise ratio (PSNR) metrics between the synCT and CT images. RESULTS: GAN training took ~11 h with a new case testing time of 5.7 ± 0.6 s. For GAN, MAEs between synCT and CT-SIM were 89.3 ± 10.3 Hounsfield units (HU) and 41.9 ± 8.6 HU across the entire FOV and tissues, respectively. However, MAE in the bone and air was, on average, ~240-255 HU. By comparison, the CNN model had an average full FOV MAE of 102.4 ± 11.1 HU. For GAN, the mean PSNR was 26.6 ± 1.2 and SSIM was 0.83 ± 0.03. GAN synCTs preserved details better than CNN, and regions of abnormal anatomy were well represented on GAN synCTs. CONCLUSIONS: We developed and validated a GAN model using a single T1-weighted MR image as the input that generates robust, high quality synCTs in seconds. Our method offers strong potential for supporting near real-time MR-only treatment planning in the brain.

Optimization of a novel large field of view distortion phantom for MR-only treatment planning.

PURPOSE: MR-only treatment planning requires images of high geometric fidelity, particularly for large fields of view (FOV). However, the availability of large FOV distortion phantoms with analysis software is currently limited. This work sought to optimize a modular distortion phantom to accommodate multiple bore configurations and implement distortion characterization in a widely implementable solution. METHOD AND MATERIALS: To determine candidate materials, 1.0 T MR and CT images were acquired of twelve urethane foam samples of various densities and strengths. Samples were precision-machined to accommodate 6 mm diameter paintballs used as landmarks. Final material candidates were selected by balancing strength, machinability, weight, and cost. Bore sizes and minimum aperture width resulting from couch position were tabulated from the literature (14 systems, 5 vendors). Bore geometry and couch position were simulated using MATLAB to generate machine-specific models to optimize the phantom build. Previously developed software for distortion characterization was modified for several magnet geometries (1.0 T, 1.5 T, 3.0 T), compared against previously published 1.0 T results, and integrated into the 3D Slicer application platform. RESULTS: All foam samples provided sufficient MR image contrast with paintball landmarks. Urethane foam (compressive strength ∼1000 psi, density ~20 lb/ft3 ) was selected for its accurate machinability and weight characteristics. For smaller bores, a phantom version with the following parameters was used: 15 foam plates, 55 × 55 × 37.5 cm3 (L×W×H), 5,082 landmarks, and weight ~30 kg. To accommodate > 70 cm wide bores, an extended build used 20 plates spanning 55 × 55 × 50 cm3 with 7,497 landmarks and weight ~44 kg. Distortion characterization software was implemented as an external module into 3D Slicer's plugin framework and results agreed with the literature. CONCLUSION: The design and implementation of a modular, extendable distortion phantom was optimized for several bore configurations. The phantom and analysis software will be available for multi-institutional collaborations and cross-validation trials to support MR-only planning.

An extended vascular model for less biased estimation of permeability parameters in DCE-T1 images.

One of the key elements in dynamic contrast enhanced (DCE) image analysis is the arterial input function (AIF). Traditionally, in DCE studies a global AIF sampled from a major artery or vein is used to estimate the vascular permeability parameters; however, not addressing dispersion and delay of the AIF at the tissue level can lead to biased estimates of these parameters. To find less biased estimates of vascular permeability parameters, a vascular model of the cerebral vascular system is proposed that considers effects of dispersion of the AIF in the vessel branches, as well as extravasation of the contrast agent (CA) to the extravascular-extracellular space. Profiles of the CA concentration were simulated for different branching levels of the vascular structure, combined with the effects of vascular leakage. To estimate the permeability parameters, the extended model was applied to these simulated signals and also to DCE-T1 (dynamic contrast enhanced T1 ) images of patients with glioblastoma multiforme tumors. The simulation study showed that, compared with the case of solving the pharmacokinetic equation with a global AIF, using the local AIF that is corrected by the vascular model can give less biased estimates of the permeability parameters (Ktrans , vp and Kb ). Applying the extended model to signals sampled from different areas of the DCE-T1 image showed that it is able to explain the CA concentration profile in both the normal areas and the tumor area, where effects of vascular leakage exist. Differences in the values of the permeability parameters estimated in these images using the local and global AIFs followed the same trend as the simulation study. These results demonstrate that the vascular model can be a useful tool for obtaining more accurate estimation of parameters in DCE studies.

A parametric model of the brain vascular system for estimation of the arterial input function (AIF) at the tissue level.

In this paper, we introduce a novel model of the brain vascular system, which is developed based on laws of fluid dynamics and vascular morphology. This model is used to address dispersion and delay of the arterial input function (AIF) at different levels of the vascular structure and to estimate the local AIF in DCE images. We developed a method based on the simplex algorithm and Akaike information criterion to estimate the likelihood of the contrast agent concentration signal sampled in DCE images belonging to different layers of the vascular tree or being a combination of different signal levels from different nodes of this structure. To evaluate this method, we tested the method on simulated local AIF signals at different levels of this structure. Even down to a signal to noise ratio of 5.5 our method was able to accurately detect the branching level of the simulated signals. When two signals with the same power level were combined, our method was able to separate the base signals of the composite AIF at the 50% threshold. We applied this method to dynamic contrast enhanced computed tomography (DCE-CT) data, and using the parameters estimated by our method we created an arrival time map of the brain. Our model corrected AIF can be used for solving the pharmacokinetic equations for more accurate estimation of vascular permeability parameters in DCE imaging studies.

Diffusion-Derived Magnetic Resonance Imaging Measures of Longitudinal Microstructural Remodeling Induced by Marrow Stromal Cell Therapy after Traumatic Brain Injury.

Using magnetic resonance imaging (MRI) and an animal model of traumatic brain injury (TBI), we investigated the capacity and sensitivity of diffusion-derived measures, fractional anisotropy (FA), and diffusion entropy, to longitudinally identify structural plasticity in the injured brain in response to the transplantation of human bone marrow stromal cells (hMSCs). Male Wistar rats (300-350g, n = 30) were subjected to controlled cortical impact TBI. At 6 h or 1 week post-injury, these rats were intravenously injected with 1 mL of saline (at 6 h or 1 week, n = 5/group) or with hMSCs in suspension (∼3 × 106 hMSCs, at 6 h or 1 week, n = 10/group). In vivo MRI measurements and sensorimotor function estimates were performed on all animals pre-injury, 1 day post-injury, and weekly for 3 weeks post-injury. Bielschowsky's silver and Luxol fast blue staining were used to reveal the axon and myelin status, respectively, with and without cell treatment after TBI. Based on image data and histological observation, regions of interest encompassing the structural alterations were made and the values of FA and entropy were monitored in these specific brain regions. Our data demonstrate that administration of hMSCs after TBI leads to enhanced white matter reorganization particularly along the boundary of contusional lesion, which can be identified by both FA and entropy. Compared with the therapy performed at 1 week post-TBI, cell intervention executed at 6 h expedites the brain remodeling process and results in an earlier functional recovery. Although FA and entropy present a similar capacity to dynamically detect the microstructural changes in the tissue regions with predominant orientation of fiber tracts, entropy exhibits a sensitivity superior to that of FA, in probing the structural alterations in the tissue areas with complex fiber patterns.

Resting state fMRI connectivity analysis as a tool for detection of abnormalities in five different cognitive networks of the brain in Multiple Sclerosis patients.

OBJECTIVES: Cognitive dysfunction is present in at least half of patients with Multiple Sclerosis. The purpose of this study was to examine functional connectivity abnormalities in patients with multiple sclerosis (MS) using resting state fMRI (rsfMRI). METHODS: Conventional MRI, rsfMRI and diffusion tensor imaging (DTI) data was acquired from 10 patients with relapsing-remitting multiple sclerosis (RRMS) and 20 healthy controls. Cross-correlation of the resting state average signal among the voxels in each brain region of the five cognitive networks: default mode network (DMN), attention, verbal memory, memory, and visuospatial working memory network, was calculated. Voxelwise analyses were used to investigate fractional anisotropy (FA) of white matter tracts. The normalized gray matter (GM), white matter and thalamus volumes were calculated. RESULTS: Compared to controls, significant deficit in MS patients at each of five networks, attention (p=0.026), DMN (p=0.004), verbal memory (p<0.001), memory (p=0.001), visuospatial working memory (p=0.003) was found. Significant reduction (p=0.034) in the normalized GM volume and asymmetry in thalamus volume (p=0.041) was detected in MS patients compared to controls. CONCLUSION: Wide spread of functional abnormalities are present within different cognitive networks in patients with RRMS, suggesting that DMN may not be sufficient for measurement of MS cognitive impairment. Larger and longitudinal studies should ascertain whether rsfMRI of cognitive networks and changes in GM and thalamus volume can be used as tools for assessment of cognition in clinical trials in MS.

Stroke increases neural stem cells and angiogenesis in the neurogenic niche of the adult mouse.

The unique cellular and vascular architecture of the adult ventricular-subventricular zone (V/SVZ) neurogenic niche plays an important role in regulating neural stem cell function. However, the in vivo identification of neural stem cells and their relationship to blood vessels within this niche in response to stroke remain largely unknown. Using whole-mount preparation of the lateral ventricle wall, we examined the architecture of neural stem cells and blood vessels in the V/SVZ of adult mouse over the course of 3 months after onset of focal cerebral ischemia. Stroke substantially increased the number of glial fibrillary acidic protein (GFAP) positive neural stem cells that are in contact with the cerebrospinal fluid (CSF) via their apical processes at the center of pinwheel structures formed by ependymal cells residing in the lateral ventricle. Long basal processes of these cells extended to blood vessels beneath the ependymal layer. Moreover, stroke increased V/SVZ endothelial cell proliferation from 2% in non-ischemic mice to 12 and 15% at 7 and 14 days after stroke, respectively. Vascular volume in the V/SVZ was augmented from 3% of the total volume prior to stroke to 6% at 90 days after stroke. Stroke-increased angiogenesis was closely associated with neuroblasts that expanded to nearly encompass the entire lateral ventricular wall in the V/SVZ. These data indicate that stroke induces long-term alterations of the neural stem cell and vascular architecture of the adult V/SVZ neurogenic niche. These post-stroke structural changes may provide insight into neural stem cell mediation of stroke-induced neurogenesis through the interaction of neural stem cells with proteins in the CSF and their sub-ependymal neurovascular interaction.

Measurement of rat brain tumor kinetics using an intravascular MR contrast agent and DCE-MRI nested model selection.

PURPOSE: Using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) in a rat glioma model, and nested model selection (NMS), to compare estimates of the pharmacokinetic parameters vp , K(trans) , and ve for two different contrast agents (CAs)-gadofosveset, which reversibly binds to human serum albumin, and gadopentetate dimeglumine, which does not. MATERIALS AND METHODS: DCE-MRI studies were performed on nine Fisher 344 rats inoculated intracerebrally with 9L gliosarcoma cells using both gadofosveset and gadopentetate. The parameters vp , K(trans) , and ve were estimated using NMS. RESULTS: K(trans) estimates using gadofosveset, compared to gadopentetate, differed in their means (gadofosveset 0.025 ± 0.008 min(-1) vs. gadopentetate 0.046 ± 0.011 min(-1) ; P = 0.0039). This difference notwithstanding, the intraclass correlation coefficient (ICC) for the two estimates of K(trans) showed nearly perfect linear dependence (ICC = 0.8479 by Pearson's r). Other estimates, ve (gadofosveset 22.7 ± 4.7% vs. gadopentetate 23.6 ± 5.6%; P = 0.4258) and vp (gadofosveset 1.5 ± 0.5% vs. gadopentetate 1.6 ± 0.4%; P = 0.25), were not different in their means between the two CAs, and there was almost perfect agreement for ve (ICC = 0.8798) and substantial agreement for vp (ICC = 0.7981) between the two CAs. CONCLUSION: Estimates of K(trans) were statistically different using gadofosveset and gadopentetate, whereas ve and vp were similar with two CAs. NMS produced robust estimates of pharmacokinetic parameters using DCE-MRI that show promise as important measures of tumor physiology and microenvironment.

Perfusion and Diffusion Abnormalities of Multiple Sclerosis Lesions and Relevance of Classified Lesions to Disease Status.

OBJECTIVE: Hemodynamic abnormality and disruption of white matter (WM) integrity are significant components in the pathophysiology of multiple sclerosis (MS) lesions. However, the roles of stratified lesions with distinct degrees of hemodynamic and structural injury in disease states remain to be explored. We tested the hypothesis that hemodynamic and structural impairment, as assessed by cerebral blood volume (CBV) and fractional anisotropy (FA), respectively, characterizes the extent of tissue injury, and the load of lesion with substantial tissue destruction would reflect the disease status and therefore, would be related to clinical disability. METHODS: Seven relapsing-remitting MS patients and seven healthy controls underwent perfusion, diffusion and conventional MRI scans. Based on T2-FLAIR and T1-weighted image, WM plaques were classified. After image coregistration, values of CBV and FA were estimated in three distinct lesion types (active, T1-hypointense and T1-isointense lesion) and compared with those obtained in WM from controls. A total of 1135 lesions were evaluated. Brain volumetric measurement and correlative analysis between brain atrophy, lesion volume and clinical disability were also performed. RESULTS: Compared with normal WM, significantly reduced CBV and FA were present in the T1-hypointense lesion, while insignificant changes in both parameters were exhibited in the T1-isointense lesion. However, increased CBV but significantly decreased FA was detected in the active lesion. A close spatial relationship between active and T1-hypointense lesion was observed. Lesion load represented by T1-hypointense plus active lesion volume significantly correlated with brain atrophy, which, in turn, significantly correlated with the severity of clinical disability. CONCLUSION: A distinct combination of CBV and FA characterizes the status of a specific lesion type. A severe structural impairment does not solely occur in the T1-hypointense lesion, but is also associated with the active lesion. The burden of the lesion with extensive structural damage provides an image index, indicative of disease status.

Characterizing brain structures and remodeling after TBI based on information content, diffusion entropy.

BACKGROUND: To overcome the limitations of conventional diffusion tensor magnetic resonance imaging resulting from the assumption of a Gaussian diffusion model for characterizing voxels containing multiple axonal orientations, Shannon's entropy was employed to evaluate white matter structure in human brain and in brain remodeling after traumatic brain injury (TBI) in a rat. METHODS: Thirteen healthy subjects were investigated using a Q-ball based DTI data sampling scheme. FA and entropy values were measured in white matter bundles, white matter fiber crossing areas, different gray matter (GM) regions and cerebrospinal fluid (CSF). Axonal densities' from the same regions of interest (ROIs) were evaluated in Bielschowsky and Luxol fast blue stained autopsy (n = 30) brain sections by light microscopy. As a case demonstration, a Wistar rat subjected to TBI and treated with bone marrow stromal cells (MSC) 1 week after TBI was employed to illustrate the superior ability of entropy over FA in detecting reorganized crossing axonal bundles as confirmed by histological analysis with Bielschowsky and Luxol fast blue staining. RESULTS: Unlike FA, entropy was less affected by axonal orientation and more affected by axonal density. A significant agreement (r = 0.91) was detected between entropy values from in vivo human brain and histologically measured axonal density from post mortum from the same brain structures. The MSC treated TBI rat demonstrated that the entropy approach is superior to FA in detecting axonal remodeling after injury. Compared with FA, entropy detected new axonal remodeling regions with crossing axons, confirmed with immunohistological staining. CONCLUSIONS: Entropy measurement is more effective in distinguishing axonal remodeling after injury, when compared with FA. Entropy is also more sensitive to axonal density than axonal orientation, and thus may provide a more accurate reflection of axonal changes that occur in neurological injury and disease.

MRI of neuronal recovery after low-dose methamphetamine treatment of traumatic brain injury in rats.

We assessed the effects of low dose methamphetamine treatment of traumatic brain injury (TBI) in rats by employing MRI, immunohistology, and neurological functional tests. Young male Wistar rats were subjected to TBI using the controlled cortical impact model. The treated rats (n = 10) received an intravenous (iv) bolus dose of 0.42 mg/kg of methamphetamine at eight hours after the TBI followed by continuous iv infusion for 24 hrs. The control rats (n = 10) received the same volume of saline using the same protocol. MRI scans, including T2-weighted imaging (T2WI) and diffusion tensor imaging (DTI), were performed one day prior to TBI, and at 1 and 3 days post TBI, and then weekly for 6 weeks. The lesion volumes of TBI damaged cerebral tissue were demarcated by elevated values in T2 maps and were histologically identified by hematoxylin and eosin (H&E) staining. The fractional anisotropy (FA) values within regions-of-interest (ROI) were measured in FA maps deduced from DTI, and were directly compared with Bielschowsky's silver and Luxol fast blue (BLFB) immunohistological staining. No therapeutic effect on lesion volumes was detected during 6 weeks after TBI. However, treatment significantly increased FA values in the recovery ROI compared with the control group at 5 and 6 weeks after TBI. Myelinated axons histologically measured using BLFB were significantly increased (p<0.001) in the treated group (25.84±1.41%) compared with the control group (17.05±2.95%). Significant correlations were detected between FA and BLFB measures in the recovery ROI (R = 0.54, p<0.02). Methamphetamine treatment significantly reduced modified neurological severity scores from 2 to 6 weeks (p<0.05) and foot-fault errors from 3 days to 6 weeks (p<0.05) after TBI. Thus, the FA data suggest that methamphetamine treatment improves white matter reorganization from 5 to 6 weeks after TBI in rats compared with saline treatment, which may contribute to the observed functional recovery.

Comparison of neurite density measured by MRI and histology after TBI.

BACKGROUND: Functional recovery after brain injury in animals is improved by marrow stromal cells (MSC) which stimulate neurite reorganization. However, MRI measurement of neurite density changes after injury has not been performed. In this study, we investigate the feasibility of MRI measurement of neurite density in an animal model of traumatic brain injury (TBI) with and without MSC treatment. METHODS: Fifteen male Wistar rats, were treated with saline (n = 6) or MSCs (n = 9) and were sacrificed at 6 weeks after controlled cortical impact (CCI). Healthy non-CCI rats (n = 5), were also employed. Ex-vivo MRI scans were performed two days after the rats were sacrificed. Multiple-shell hybrid diffusion imaging encoding scheme and spherical harmonic expansion of a two-compartment water diffusion displacement model were used to extract neurite related parameters. Bielshowski and Luxol Fast blue was used for staining axons and myelin, respectively. Modified Morris water maze and neurological severity score (mNSS) test were performed for functional evaluation. The treatment effects, the correlations between neurite densities measured by MRI and histology, and the correlations between MRI and functional variables were calculated by repeated measures analysis of variance, the regression correlation analysis tests, and spearman correlation coefficients. RESULTS: Neurite densities exhibited a significant correlation (R(2)>0.80, p<1E-20) between MRI and immuno-histochemistry measurements with 95% lower bound of the intra-correlation coefficient (ICC) as 0.86. The conventional fractional anisotropy (FA) correlated moderately with histological neurite density (R(2) = 0.59, P<1E-5) with 95% lower bound of ICC as 0.76. MRI data revealed increased neurite reorganization with MSC treatment compared with saline treatment, confirmed by histological data from the same animals. mNSS were significantly correlated with MRI neurite density in the hippocampus region. CONCLUSIONS: The present studies demonstrated that neurite density can be estimated by MRI after TBI and MRI measurement of neurite density is a sensitive marker to MSC treatment response.

MRI detects brain reorganization after human umbilical tissue-derived cells (hUTC) treatment of stroke in rat.

Human umbilical tissue-derived cells (hUTC) represent an attractive cell source and a potential technology for neurorestoration and improvement of functional outcomes following stroke. Male Wistar rats were subjected to a transient middle cerebral artery occlusion (tMCAo) and were intravenously administered hUTC (N = 11) or vehicle (N = 10) 48 hrs after stroke. White matter and vascular reorganization was monitored over a 12-week period using MRI and histopathology. MRI results were correlated with neurological functional and histology outcomes to demonstrate that MRI can be a useful tool to measure structural recovery after stroke. MRI revealed a significant reduction in the ventricular volume expansion and improvement in cerebral blood flow (CBF) in the hUTC treated group compared to vehicle treated group. Treatment with hUTC resulted in histological and functional improvements as evidenced by enhanced expression of vWF and synaptophysin, and improved outcomes on behavioral tests. Significant correlations were detected between MRI ventricular volumes and histological lesion volume as well as number of apoptotic cells. A positive correlation was also observed between MRI CBF or cerebral blood volume (CBV) and histological synaptic density. Neurological functional tests were also significantly correlated with MRI ventricular volume and CBV. Our data demonstrated that MRI measurements can detect the effect of hUTC therapy on the brain reorganization and exhibited positive correlation with histological measurements of brain structural changes and functional behavioral tests after stroke. MRI ventricular volumes provided the most sensitive index in monitoring brain remodeling and treatment effects and highly correlated with histological and functional measurements.

MRI measurement of angiogenesis and the therapeutic effect of acute marrow stromal cell administration on traumatic brain injury.

Using magnetic resonance imaging (MRI), the present study was undertaken to investigate the therapeutic effect of acute administration of human bone marrow stromal cells (hMSCs) on traumatic brain injury (TBI) and to measure the temporal profile of angiogenesis after the injury with or without cell intervention. Male Wistar rats (300 to 350 g, n=18) subjected to controlled cortical impact TBI were intravenously injected with 1 mL of saline (n=9) or hMSCs in suspension (n=9, 3 × 10(6) hMSCs) 6 hours after TBI. In-vivo MRI acquisitions of T2-weighted imaging, cerebral blood flow (CBF), three-dimensional (3D) gradient echo imaging, and blood-to-brain transfer constant (Ki) of contrast agent were performed on all animals 2 days after injury and weekly for 6 weeks. Sensorimotor function and spatial learning were evaluated. Volumetric changes in the trauma-induced brain lesion and the lateral ventricles were tracked and quantified using T2 maps, and hemodynamic alteration and blood-brain barrier permeability were monitored by CBF and Ki, respectively. Our data show that transplantation of hMSCs 6 hours after TBI leads to reduced cerebral atrophy, early and enhanced cerebral tissue perfusion and improved functional outcome compared with controls. The hMSC treatment increases angiogenesis in the injured brain, which may promote neurologic recovery after TBI.